The two Josephson junction flux qubit with large tunnehng amphtude 
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In this paper we discuss solid-state nanoelectronic realizations of Josephson flux qubits with large 
tunneling amplitude between the two wells. The latter can be controlled via the height and form of 
the barrier, which is determined by quantum-state engineering of the flux qubit circuit. The simplest 
circuit of the flux qubit is a superconducting loop interrupted by a Josephson nanoscale tunnel 
junction. The tunneling amplitude between two macroscopically different states can be essentially 
increased, by engineering of the qubit circuit, if tunnel junction is replaced by a ScS contact. 
However, only Josephson tunnel junctions are particularly suitable for large-scale integration circuits 
and quantum detectors with preset-day technology. To overcome this difficulty we consider here 
novel flux qubit with high-level energy separation between "ground" and "excited" states, which 
consists of a superconducting loop with two low-capacitance Josephson tunnel junctions in series. We 
demonstrate that for real parameters of resonance superposition the splitting amplitude can reach 
values greater than 1 K. Analytical results obtained with instanton technique show good correlation 
with a numerical exact solution. 

PACS numbers: 03.75.Lm, 74.50.+r, 85.25.Cp 



Since successful demonstration of Rabi oscillations and 
Landau-Zener coherent effects,— '^'^'^'^ the superconduct- 
ing qubits (quantum bits) based on mesoscopic Josephson 
junctions became the subject of consideration as possible 
candidates to be the basic elements of a quantum com- 
puter hardwarCj^'" including detectors to measure the 
state of an individual qubit i^i^ i^^d^ The Josephson junc- 
tion (JJ) qubits have two energy scales which are the 
Josephson coupling energy Ej and the charging energy 
Ec of the JJ, and they are subdivided into flux qubits, 
charge qubits, as well as charge-phase and persistent- 
current qubits. In principle, all circuits of a quantum 
computer can be fabricated by modern techniques using 
these superconducting qubits. However, it is but poor 
qualitySii^. of the experimentally tested elements that is 
the limiting factor on the way of implementation of quan- 
tum registers. For example, an important but still un- 
solved problem in the physics of a qubit working in the 
charge regime with Ec/Ej 3> 1 is an essential decrease 
of high spectral density of the noise associated with the 
motion of charge in traps. In its turn, the phase qubit 
{Ej/Ec ^1), which utilizes the phase of the supercon- 
ducting order parameter as a dynamic variable, is much 
less sensitive to the charge fluctuations but is subject to 
the influence of the noise in critical current of J J, spin 
fluctuations and Nyquist noise currents generated by ex- 
cess ambient temperature. The tunnel splitting of the 
energy levels arising from the coherent superposition of 
the macroscopic states is small usually, AEqi ~ 150 - 250 
mK. Taking into account the effective noise temperature, 
which can reach Te// ~ 50 - 100 mK in experimental 
studies of the qubit dynamics, leads to -'^^d'^d^ a dramatic 
fall of the decoherence times and relaxation times . 
This means that, in order to enhance considerably the 



qubit quality^, (the number of one-bit operations during 
the coherence timespan), the system with large {AEqi 
> ^O^^'^ - 10~^^ J) tunnel splitting of the energy levels 
should be created. 

Undoubtedly, the problem of creation of a quantum 
register based on Josephson qubits brings up many is- 
sues but presently the invention of a high-quality qubit 
is the most important one among them. It is easy to show 
that the rate of the energy exchange between two macro- 
scopic states in a flux qubit is bounded by the "cosine" 
shape of the potential barrier and cannot be increased 
owing to decreasing the barrier height since the latter de- 
termines the characteristic rate of thermal decay of the 
current-flow states. A similar limitation associated with 
the lowering of the effective barrier height can appear also 
when highly increasing the pre-exponential factor. It is 
absolutely obvious that the ideal case for a flux qubit is 
when the tunnel barrier in the phase space looks like 11- 
function having sufficiently large height and small action. 
It was this issue that motivated the authors of the Ref. 
[T^ for analyzing the phase-slip qubit, whose creation re- 
quired developing a new non-Josephson technology. In 
this paper we search for an improved barrier design for 
the J J flux qubit. 

The recent Ref. [13 demonstrated how the level split- 
ting can be increased at low temperatures (T — > 0) by 
an order of magnitude with the potential barrier height 
kept unchanged by modifying the qubit 's potential bar- 
rier shape due to using the clean-limit ScS junction in the 
superconducting ring. However, the fabrication difficul- 
ties of obtaining pure and reproducible ScS junctions are 
the serious hindrance in the way of designing large-scale 
integrated qubit circuits. 

To solve this problem, the analysis is carried out in the 
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FIG. 1: Schematic picture of the proposed 2JJ flux qubit 
with a SQUID configuration (a) and its circuit diagram (b). 
The loop carrying supercurrent Is is pierced by an externally 
applied magnetic flux <[>e (towards an observer). The indi- 
vidual SIS Josephson junctions are characterized by coupling 
energies E,j\,Ej2, critical currents Ici,Ic2 and capacitances 
Ci , C2 which do not difi'er significantly. The loop induc- 
tance L is small enough so that the 2JJ SQUID has only two 
metastable flux states. The parameter g = Eji/Eci ^ 1. 



present paper of the two-Josephson-junction flux qubit 
(2 J J flux qubit), which can be considered as a super- 
conducting ring of inductance L interrupted by two al- 
most equivalent tunnel SIS junctions with the Joseph- 
son energies Eji,Ej2, the critical currents Ici,Ic2 and 
the capacitances Ci,C2, respectively (see Fig. [I{a,b)). 
The difference between the two SIS mesoscopic junc- 
tions will be characterized by the asymmetry parame- 
ter A = Ic2/Ici = Ej2/Eji = C2/C1 < 1 so that the 
"junction 1" would have greater or equal values of the 
Josephson energy, the critical current and the capaci- 
tance as compared to the "junction 2" . The external 
magnetic flux <i>e can be coupled to the qubit by a sep- 
arate coil located in close proximity to the qubit's loop. 
It is well known that in classical limit the circulating 
current Is as a function of external magnetic flux for dc 
SQUIDs with Id = Ic2 has the singularity in the points 
of <I>e = '^oin + 1/2) ($0 is the flux quantum) so that 
the two-junction interferometer can be considered as a 
"single-junction" one, with the potential energy shape 
in the phase space being modified. To retain this classi- 
cal effect in the quantum regime the proposed 2JJ flux 
qubit have to satisfy the usual condition for flux qubits 
g = Eji/Eci 3> 1, at which phase is a good quantum 
variable and the charging effect on the island between JJ 
contacts is negligiblei^ 

The problem lies in determining and analyzing the tun- 
nel splitting AE'oi — Ei — Eq of the degenerated zero 
energy level in the double-well symmetric potential of a 
2JJ flux qubit (at corresponding external conditions) re- 
sulted from the coherent quantum tunneling of the mag- 
netic flux between the wells. In the proposed mesoscopic 
system in quantum regime, the two lower energy levels Eq 
and El arising from coherent superposition of the macro- 
scopically distinct flux or persistent-current states form 
a quantum bit (qubit). It turns out that, because of the 



change in the form of the potential energy of the 2JJ flux 
qubit as compared to the IJJ qubit, the tunnel splitting 
AEqi can multiply rise reaching the values > 1 (in tem- 
perature units) and substantially enhance the properties 
of the qubit as a basic element for quantum computa- 
tions. However, even slightly away from the symmetry 
condition Id = Ic2, the splitting energy is reduced. The 
sensitivity of the AEqi magnitude to A as well as to the 
junction parameters can limit applications based on the 
2JJ flux qubit both for quantum computation and quan- 
tum detectors. 

With this sake, we will discuss the 2JJ flux qubit in 
the approximation of the Hamiltonian of an isolated sys- 
tem in the zero temperature limit. All the dissipative 
processes associated with the own and the external, re- 
garding the system, degrees of freedom (the quasiparti- 
cles, the magnetic flux fluctuations in the qubit and in 
the outer measuring circuit, etc.) arc neglected in this 
approximation. In the framework of this approximation, 
only the supercurrent component flows in the qubit ring 
which in classical regime, according to the Josephson re- 
lation, is equal to 

Is = Id sin (fi = Ic2 sin (^2 , (1) 

where tpi , (^2 are the order parameter phase differences at 
corresponding tunnel junctions. It is convenient to count 
the values of the supercurrent Is and the phase differences 
at the junctions clockwise, the applied magnetic flux $e, 
the total magnetic flux in the ring and the supercurrent 
Is being tied by the relation $ = <I>e — L/s($) . The classic 
Hamiltonian of the 2JJ flux qubit in the approximation 
of the isolated system contains the contributions of the 
electrostatic energy of the charges in the junction capaci- 
tances, the junction Josephson energies and the magnetic 
energy of the supercurrent in the ring, and has the form: 

- {Eji cos ipi + Ej2 cos (P2) + — h Eq , 

where A^o is the number of the excess (deficient) Cooper 
pairs in the banks of the SIS Josephson junctions, Eq 
is the constant fixing the reference level for the po- 
tential energy. Using relation we will reduce ex- 
pression for the Josephson energy in Hamiltonian ^ 
to the form Ej{(f)) = —{Ejicosipi + Ej2COSip2) — 
—Ej\\/ (1 — A)2 -I- 4Acos^(0/2) , where a new variable of 
the overall phase = Lp\ ~\- Lp2 is introduced. As follows 
from the theory of charge-phase qubits*'^^ this classical 
expression Ej((f>) for the Josephson energy of the two- 
junction interferometer remains unchanged in the quan- 
tum regime of system dynamics in the limit g ^ 00. 
Mathematically it means a fulfilment of the relationship 
(fT|) between semiclassical variables ifi , ip2 in the quantum 
regime at g ^ 1, and physically disappearance of phase 
fluctuations induced by the charging effect on the island 
between JJ contacts^ 
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Due to the unambiguity of the superconductor's order 
parameter the variable cj) satisfies the condition 



= 27r 

$0 



27r7i, $ = n$oH $0, ^0 = 7rS/e , 

27r 

(3) 

where n is the integer number of the flux quanta $0 in 
the total magnetic flux $ (below, we will consider the 
qubit to work in the n — mode). Owing to the re- 
lationships H]), ([3]) there is the only independent phase 
variable from ipi,ip2,(t>, and in the quantum regime the 
physical fluctuating variable is the total phase difference 
(j) at the both junctions which equals, to within 2n fac- 
tor, to the total magnetic flux in the ring in the units of 
flux quantum (0/27r = $/$o)- The other two variables, 
and </52 are unambiguously depend on (f>, i.e. 
the total phase difference (j) uniquely expands into the 
components ipi , ip2 (see below) . 

The transition to the quantum description of the 
flux qubit consists in associating the value A^o of the 
Cooper pairs tunneling through the junctions with the 
d 

operator A'o = —i-^ , conjugated to the phase oper- 
dcp 

ator (j) ([No,(j)\ — —i), and solving the Schrodinger 
equation with the obtained quantum Hamiltonian in 0- 
representationjii By applying the quantization procedure 
to the Hamiltonian Q and writing down the energy con- 
tributions via the variable 0, we will come to a canoni- 
cal form of the Hamiltonian of the 2JJ flux qubit in the 
quantum case: 

^=2M+^(^)--2Ma^+ 



+E. 



71 £0 



(l-A)2+4Acos2. 



2Pl 



(4) 

which can be considered as the Hamiltonian of a quantum 
particle with the mass M moving in the potential U{4>). 
d . - 

Here P = HNq — ~ifi— ([P, 0] = —ih) corresponds to 
aq) 

_ f2TTV ACi 

the particle momentum operator, M — — — 

V^'o/ (A + 1) 

27r /27rV 
is its mass, Pl — — Lid — — I LEji is the po- 

tential parameter, 0e = 2n — is the external magnetic 

flux parameter (the constant Eq is chosen further from 
the condition for the symmetric potential of equalling to 
zero in the minima points). The flux qubit parameter 

p2 



2Ci Ci 



g — Ejl/Eci r 1 o , 

The potential J7(0) shape depends on the parameters 
A, PL,4'e- We are interested in the case of symmetric po- 
tential which, according to (|4]), is realized at 0e = tt 
(<I>e = ^"0/2). It should be noticed that in the extreme 
case of identical junctions, at A = 1, the potential C/(0) 
coincides with the potential of the flux qubit based on 
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FIG. 2: Potential U{i^/^o)/kB in temperature units for 
IJJ qubit with /3l = 1.602 - 1 and for 2JJ qubit for sev- 
eral parameter couples (A,/3l): (0.99,0.9) - 2, (0.9,1.058) 
- 3, (0.8,1.276) - at external magnetic flux <I>e = $o/2. 
The potential barrier heights Uo in curves 1-4 a-re equal, 
Uo/kB = 9.64 . 



the clean ScS contact studied in Ref. 10 (Equation (3)): 



UscsW^Ejl -2 \cos{^/2)\ 



2(3l 

At the same time, owing to renormalization of the mass 
M for the 2JJ flux qubit by the factor A/(A + 1) in re- 
spect to the corresponding mass for the ScS flux qubit 
(provided that the capacitances of SIS and ScS junctions 
are equal, Ci = C), for A = 1, the relation for masses 
is M2J.J = Mscs/"^ ■ Hence, at A w 1 the splitting Ai?oi 
in 2JJ qubit is expected yet more than in the ScS qubit. 
The parameter (3^ determines the height of the potential 
barrier of the double-well potential, so that the barrier 
height goes down while reducing (3^. Like in the case of 
ScS qubit, the 2JJ qubit potential has two local minima 
even at < 1 (unlike the SIS qubit where the double-well 
potential exists at /? > 1 only), which gives a possibility 
of considerable scaling down the geometric dimension (in- 
ductance) of the system with the mesoscopic junctions. 

Fig. [2] shows the potential f7($/$o)/fci3 of 2JJ 
flux qubit for several parameters couples (A,/3l) and 
also, for the comparison sake, the well-known potential 
Usis{'^/'^o)lkB of IJJ flux qubit at external magnetic 
flux $e — $0/2. The inductances L for both types of 
the qubits can be supposed equal (to specify the mag- 
netic flux fluctuation level) while the parameter [3l (i.e. 
the critical currents /ci, Ic of the corresponding SIS junc- 
tions) in all the dependences are chosen in such a way so 
that the potential barriers in all the potentials were of the 
same height [/q- The latter requirement implies roughly 
equal decay rates for the metastable states due to thermal 
fluctuations, taking them into account being beyond our 
consideration. Apparently, to realize the quantum regime 
in a physical experiment, the value Uo/ks must highly 
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exceed the system temperature. As seen from Fig. [2l the 
potentials U{^/^o)/kB for a 2JJ qubit have lesser width 
(between the potential minima points) as compared to 
the corresponding potentials for a IJJ qubit while the 
area under the potential curve between the points of its 
minima for the 2JJ qubit shrinks greatly against the cor- 
responding area for the IJJ qubit. Additionally, if the 



corresponding capacitances of the SIS junctions in both 
2JJ and IJJ qubits are equal (Ci = C) then the ratio 
of the effective masses for these qubits is A/(A + 1). As 
it will be shown below, it is the change in the potential 
shape and the decrease of the effective mass in 2JJ qubit 
that lead to multiple rise in the amplitude of its tunnel 
splitting. 




FIG. 3: a) Integral phase-current relation Is{(p/2Ti)/ Ic\ for 
2JJ qubit at various A; 0.99 - 1, 0.9 - 2, 0.8 - 3, 0.5 - I 
b) Functions [ipi /2-K)[(f>/2-n:] - 1, (<^2/27r)[(/>/27r] - 2 for 2JJ 
qubit at A = 0.9. The straight line 1/^1(0) -|- V'2(</>) = <j> - 
3 corresponds to the (f) definition. The values of ip\/2n — 
arcsin(A)/27r « 0.18 and ip2/2n — 0.25 (the latter being 
A-independent) correspond to (f)m/2n — arccos(— 0.9)/27r ~ 
0.43. 



The shape of the potential ?7(0) for 2JJ qubit is di- 
rectly relates to the form of the current-phase relation 
derived from H]), dS]): 



Id 



~ sin ipi — X sin Lp2 = 



A sin ( 



v/(l - A)2 +4Acos2(0/2) 



(5) 

The current-phase relation Is{4>) extrema (which are 
equal by their absolute values) are located in the points 
(f)rn = arccos(— A) (maximum; 7r/2 < </),„ < vr) and 
(j}mi = 27r— arccos(— A) (minimum) symmetrically around 
the point </) = tt, where the supercurrent vanishes to 
zero (Js ~ 0) alternating its direction. Thus, at A be- 
ing near the unity, in the interval {(j),m (j^mi), the super- 
current Is changes dramatically from its maximum to 
minimum value with alternating the current direction in 
the point cj) = n. Fig. ^a.) displays the integral current- 
phase dependence /s(0/27r)//ci for 2JJ qubit for several 
parameters A. The interval {4'm,4'mi) shrinks as the pa- 
rameter A increases and the maximum-to-minimum by- 



current transition becomes more sharp. The extreme 
case A = 1 corresponds to — (f>mi — tt with the infi- 
nite derivative of the current-phase relation in the point 
TT. Let us also consider the order parameter phase dif- 
ferences (^1(0), ip2{4'), derived directly from The 
analysis of formula ^ shows that the function <^i(<^) 
for a junction with high critical current has extrema in 
the points 'prmi'mi- The transition from the maximum 
positive value 'pi{(t>m) — arcsinA (0 < fi{4>m 
to the minimum negative value yi((/'mi) = —arcsinA 
with alternating the phase difference sign in the point 
TT (ifi^Tr) = 0) takes place in the interval {4>mi (t>mi), and 
fiiO) = (/?i(27r) = 0. The function (/32(0) for a junction 
with lower critical current is a monotonically increasing 
one from <y32(0) = to (^2(27r) = 2tt, which is symmet- 
rical with respect to the line y — (j)] '^2{'^) = tt, and 
f2{4'm) = 7r/2, ip2i4'mi) = 37r/2. In the extreme case of 
A = 1 the functions ipi{4>), ^2(4') behave as follows: ipi — 
f2 ~ 0/2 at < (/) < tt; at the point tt a jump appears 
in the function VJi(0) between the values 7r/2, — 7r/2 with 
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further linear rise up to (Pi{2tt) = 0, while the function 
f2i4') demonstrates a jump between the values 7r/2, 37r/2 
with further linear increase up to ip2{27T) = 27r. Fig. El^b) 
exhibits dependences ((^i/27r)[0/27r], {lp2/27t)[(I)/2tt] for a 
certain A, with their distinctive appearance. The straight 
line ipi{(j)) + ^2(4') — 4' corresponds to the <f> definition 
showing the expansion of the total phase difference over 
the both junctions into the component phase differences 
of the order parameter over each of them. 

We will find the tunnel splitting AiJgi of the degener- 
ated zero level in the symmetrical (at (pe — tt) double- 
well potential in the 2JJ flux qubit by numerical 
solution of the Schrodinger equation and analytically by 
using instanton technique in the semiclassical approxi- 
mation. To find a numeric solution of the Schrodinger 
stationary equation 



(6) 



with Hamiltonian Q, a kind of the finite elements 
method is used with approximation of the potential ?7(0) 
by a piecewise constant function. Zero boundary condi- 
tions are used for the wave function ^'(0), the domain 
width and the element quantity being set so that provide 
good accuracy of the calculation. 

In the semiclassical approximation the problem of a 
tunneling quantum particle can be solved using the in- 
stanton technique J^ii^ For a particle of the mass Af, mov- 
ing at zero temperature in symmetric double-well poten- 
tial Vix), referenced from its minimum level (y(±a) = 0, 



±a are the minimum points), the expressions for the en- 
ergy levels E\_Q and the tunnel splitting Ai?oi read like 



£^1,0 = £^0 ± — ^ = — ± fiifexp I - Y 



^0 = / dx^J2MV(x) 



(7) 



M 



AEqi = 2A hujQ \l exp I — — 



Here luq is the frequency of the particle zero oscillations 
in each of the wells, 5*0 is the particle action on the in- 
stanton trajectory, the dimensionless constant A is found 
from the equation for the instanton's function t{x) in the 
asymptotic limit: 



^('^) la; — ta 



dx 



1 , a — X 
= in ■ 



y/2V{x)/M ojo Aa ' 



A=limi^exp( r..,/^^^-o 



(8) 



2V{x) 



Starting from the Hamiltonian Q and using formulae 
(O,®, we obtain the tunnel splitting Ai?oi for the 2JJ 
flux qubit in the case of symmetric double- well potential: 




A variable a = — 7r)/2 is introduced in formulae 
(O (due to the potential symmetry condition 0e = tt), 
the minima point of ao > of the potential Vipl) sat- 
isfying equation ((9^). The accuracy of the semiclassi- 
cal approximation is high provided that S'o/fi ^ 1, the 



method accuracy degrades as the dimensionless variable 
S^jh diminishes approaching the unity. The results of a 
numerical analysis is of great importance in this region. 
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FIG. 4: a) Function AEoi{l3L)/kB for 2JJ qubit at A = 0.99 
- 1, ScS qubit - 2, IJJ qubit - 3; the points in the numerical 
curves corresponding to equal height (9.64 K) of the poten- 
tial barrier for all the qubits are indicated by arrows, b) 
Function Ai5oi(/3L)/A:B for 2JJ qubit and various A: 0.99 - 
1, 0.9 - 2, 0.8 - 3 and "level line" of equal heights (9.64 K) 
of the potential barriers at varying X - 4- The numerically 
obtained results are pointed by hollow circles, the analytically 
obtained ones are plotted by solid lines in the graphs (a) and 
(b). Dashed lines show the lowest boundary Pl, at which the 
level El height becomes equal to the potential barrier height 
Uo- For IJJ and ScS qubits, the capacitance of corresponding 
(SIS and ScS) junctions is C = 2.7-10^^^ F while for 2JJ qubit 
the capacitance of the larger SIS junction is Ci = 2.7T0~^^ F. 
The geometric ring inductance is L = 3.0-10"^'' H, parameter 
g « 76/3l for all the qubits. 



Fig. |ll^a,b) presents the /3L-dependences of the tunnel 
splitting AEoi{f]L)/kB for 2JJ, ScS and SIS flux qubits 
at the equal capacitances of the corresponding junctions 
Ci = C = 2.7 fF and at the inductance L = 0.3 nH of the 
qubits loop. In both plots the curves calculated numeri- 
cally are pointed by hollow circles while the ones obtained 
analytically using the instanton technique are plotted by 
solid lines. The formulae (O were used for 2JJ qubit 
while similar formulae were taken for ScS and SIS qubits 
based on the forms of their potentials. The change in the 
parameter Pl means the variation of the critical currents 
Ici,Ic of the corresponding junctions at a fixed induc- 
tance L. The double- well potential height decreases with 
lowering the parameter Pl, the energy level Ei being 
equalized to the potential barrier height Uq at a certain 
Plo {El = Uo) and exceeding it with further (3l lower- 
ing. Then, the wave function corresponding to the level 
El is no further a superposition of the states localized in 
the left and right wells. The boundary values (3lo ior the 
curves in the figure are indicated by dash lines. In the 
vicinity of Plo, a-t (Uq — Ei) ksT , the quantum coher- 
ence will be destroyed due to thermal fluctuation causing 
the over-barrier transitions. One can see from Fig. |3]that 



the numerically and the analytically obtained curves al- 
most coincide at large Pl and begin to diverge at lower 
Pl- This is because of the condition of semiclassicity 
So/h^ 1 starts to fail when diminishing P^. This, in its 
turn, is caused by decreasing of the barrier height Uq and 
therefore the action Sq. The analysis of the dependences 
So{PL)/fi reveals that So/h ^ 1 at /3 ~ Plo and the 
relative divergence between the numerical and the ana- 
lytical results for 2JJ and ScS qubits is within 2 to 10 per 
cent. For a SIS qubit a flt of the numerical and analytical 
results requires the more accurate fulfilment of the semi- 
classicity condition. However, it is follows even from this 
analysis that obtaining the tunnel splitting AEqi > 1 K 
in the flux qubit based on a single SIS junction is impos- 
sible under condition of weak {Uq — i?i ^ ksT) influence 
of thermal fluctuations on the metastable states decay. 
The dependences So{PL)/fi for 2JJ, ScS and SIS qubits 
are close to linear ones, whose slope (the action Sq from 
Pl rate of increase) being higher in the indicated order. 
The value of tunnel splitting in the region of its expo- 
nential smallness 5*0 (/3l)/^ ^ 1 diminishes in the same 
sequence. 

The points in the numerical curves corresponding to 
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the equal heights of the potential barriers {Uq = 9.64 K) 
are indicated by arrows in Fig. [31^ a). The correspond- 
ing values of the parameter couples £^oi(/3L)/fcB) for 
2JJ (A = 0.99), ScS and SIS flux qubits are: (0.9, 3.83 K), 
(0.88, 1.79 K), (1.60, 0.16 K). It is seen that, under this 
condition, the tunnel splitting in a 2 JJ qubit is more than 
twice the splitting in a ScS qubit and more than 20 times 
higher than that of a SIS qubit. The curve for the tun- 
nel splitting in 2JJ lies completely above the curves for 
ScS and SIS qubits, and the tunnel splitting for a 2JJ 
qubit reaches the value of3.14K at Pl = 1 > /3io- The 
advantages of a ScS qubit if compared to a SIS qubit 
was thoroughly analyzed in Ref. [lO|. Note that yet more 
increase of the tunnel splitting in a 2JJ qubit in com- 
parison with a ScS qubit with the matched parameters 
mentioned above results from the fact that their poten- 
tials (at A w 1) practically coincide while the effective 
mass M in 2JJ qubit is less by a factor of about two. 
Fig. m^b) shows the dependence Ai?oi(/3L)/fcB for the 
2JJ qubit at several A, and also a "level line", the line 
A£^oi(/3i)/fcB corresponding to equal height (9.64 K) of 
the potential barriers in 2JJ qubit with varying A. The 
curve AEoi{PL)/kB shifts right when decreasing A, and 
the smaller being the value of A, the higher the tunnel 
splitting at a fixed (3l- This, however, is due to the low- 
ering of the barrier Uq height when decreasing A that 



leads to the exponential rise of the thermal decay rate. 
Note that when desymmetrizing the junctions a fit be- 
tween the numerical and the analytical curves gets worse 
because Sq{I3l)/^ decreases. As seen from the plot, the 
value of the tunnel splitting gradually diminishes while 
moving along the level line with the equal height of the 
potential barriers towards the lower values of the junction 
symmetry parameter A (and the higher [3l)- 

It should be emphasized that the principal require- 
ments to 2JJ flux qubits, namely: A > 0.95; C < 
SOfF/^m^; jc ~ 10^ A/cm^; 1^ ^ l^iA at the JJ area 
Sj ~ 0.1 IMS?] L ~ 0.3 nH, /J^ ~ 1 can be met with the 
present-day technology based on Nb, NbN, MoRe mate- 
rials with superconductivity gap A(0) ~ lOK (see, e.g., 
Ref. '20'). One can notice in conclusion that 2JJ flux 
qubit with large amplitude of tunnel splitting potentially 
has some strong advantages: (i) weak sensitivity to the 
motion of charge in traps; (ii) extremely fast excitation 
(pumping frequency) in qubit-based readout as well as in 
computer circuits due to considerable increasing of the 
quantum tunneling rate v ~ AiSgi; (iii) macroscopically 
large energy relaxation times (see, e.g., Ref. and 
Refs. therein); (iv) further improvement of qubit coher- 
ence characteristics;^ 

We thank W. Krech and D. Born for helpful discus- 
sions. 



* Electronic address: shnyrkov@ilt.kharkov.ua 

^ Y. Nakamura, Yu.A. Pashkin, and J.-S. Tsai, Nature 398, 
786 (1999). 

^ D. Vion, A. Aassime, A. Cottet, P. Jooyez, H. Pothier, C. 
Urbina, D. Esteve, and M.H. Devoret, Science 296, 886 
(2002). 

3 Y. Yu, S. Han, Xi Chu, S.-I. Chu, and Z. Wang, Science 
296, 889 (2002). 

* J.M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 
Phys. Rev. Lett. 89, 117901 (2002). 

^ M. SiUanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and P. 
Hakonen, cond-mat/0510559 vl (2005). 

^ G. Wendin and V. S. Shumeiko, "Supercomducting cir- 
cuits, qubits and computing" , in Handbook of Theoretical 
and Computational Nanotechnology, (American Scientific 
Publishers, ed. M. Rieth and W. Schommers, 2006). 
Y. Makhlin, G. Schon, A. Shnirman, Physica C 368, 276 
(2002). 

* A.B. Zorin, Phys. Rev. Lett. 76, 4408 (1996); A.B. Zorin, 
Physica C 368, 284 (2002). 

^ V.I. Shnyrkov, G.M. Tsoi, D.A. Konotop, and I.M. 
Dmitrenko, in Proc. 4th Int. Conf. SQUID'91 (Session on 
SET and Mesoscopic Devices), p. 211, Berlin, Germany 
(1991). 

^° V.I. Shnyrkov, S.I. Melnik, Low Temp. Phys. 33, v. 1, 15 



(2007). 

" E. n'ichev, Ya. S. Greenberg, Europhys. Lett. 77, 58005 



( 2007) ; |cond-mat /06084 1 6 

J.E. Mooij and C.J. P.M. Harmans, New J. Phys. 7, 219 
(2005). 

" M. Grifoni, E. Paladino, and U. Weiss, Eur. Phys. J. B 10, 
719 (1999). 

^'^ L. Tian, S. Lloyd, and Orlando, Phys. Rev. B 65, 144516 
(2002). 

A.Yu. Smirnov, Phys Rev. B 67, 155104 (2003); Phys Rev. 
B 68, 134514 (2003). 

W. Krech, M. Grajcar, D. Born, I. Zhilyaev, Th. Wagner, 
E. Il'chev, and Ya. Greenberg, Phys. Lett. A 303, 352 
(2002). 

A.J. Leggett, J. Phys.: Condens. Matter 14, R415 - R451 
(2002). 

S. Coleman, Aspects of symmetry - selected Erice lectures 
(CUP 1985), Ch. 7. 

H. Kleinert, Path Integrals in Quantum Mechanics, Statis- 
tics, and Polymer Physics, and Financial Markets (World 
Scientific Publishing Company, 3 edition, 2004), Ch. 17. 
H. Zangerle, J. Konemann, B. Mackrodt, R. Dolata, S.V. 
Lotkhov, S. A. Bogoslovsky, M. Gotz, and A.B. Zorin, 
Phys. Rev. B 73, 224527 (2006). 



